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Abstract. Using a generalization of the Madelung transformation, we derive the 
hydrodynamic representation of the Klein-Gordon-Einstein equations in the weak field limit. We 
consider a complex self-interacting scalar field with an arbitrary potential of the form U(|</?| 2 ). 
We compare the results with simplified models in which the gravitational potential is introduced 
by hand in the Klein-Gordon equation, and assumed to satisfy a (generalized) Poisson equation. 
Nonrelativistic hydrodynamic equations based on the Schrodinger-Poisson equations or on the 
Gross-Pitaevskii-Poisson equations are recovered in the limit c —¥ + 00 . 


1. Introduction 

Scalar fields (SF) play an important role in particle physics, astrophysics, and cosmology [T]. 
Their evolution is usually described by the Klein-Gordon (KG) equation [2] [3] which can be 
viewed as a relativistic extension of the Schrodinger equation [3J. The KG equation describes 
spin-0 particles (bosons) that can be charged (complex SF) or neutral (real SF). The coupling 
between the KG equation and gravity through the Einstein equations, leading to the Klein- 
Gordon-Einstein (KGE) equations, was first considered in the context of boson stars [5, 3 E]- 
It has also been proposed that dark matter (DM) halos may be made of a SF described by 
the KGE equations (see, e.g., isua m for recent reviews). Actually, at the galactic scale, 
the Newtonian limit is valid so that DM halos can be described by the Schrodinger-Poisson 
(SP) equations or by the Gross-Pitaevskii-Poisson (GPP) equations. In that case, the wave 
function t ) describes a Bose-Einstein condensate (BEC) at T = 0, and the self-interaction 

of the bosons is measured by their scattering length a s . Therefore, DM halos could be gigantic 
quantum objects made of BECs. The wave properties of bosonic DM may stabilize the system 
against gravitational collapse, providing halo cores and sharply suppressing small-scale linear 
power. This may solve the problems of the cold dark matter (CDM) model such as the cusp 
problem and the missing satellite problem. The scalar field dark matter (SFDM) model and the 
BEC dark matter (BECDM) model, also called $DM models, have received much attention in 
the last years. 

Since DM may be a SF, it is of considerable interest to study the cosmological implications 
of this scenario. The cosmological evolution of a spatially homogeneous noninteracting real SF 
described by the KGE equations competing with baryonic matter, radiation and dark energy 


was considered by Matos et al. m ■ They found that real SFs display fast oscillations but 
that, on the mean, they reproduce the cosmological predictions of the standard ACDM model. 
The study of perturbations was considered by Suarez and Matos [12] for a self-interacting 
real SF described by the Klein-Gordon-Poisson (KGP) equations and by Magana et al. [13] 
for a noninteracting real SF described by the KGE equations. These studies show that the 
perturbations can grow in the linear regime, leading, in the nonlinear regime, to the formation 
of structures corresponding to DM halos. This is in agreement with the early work of Khlopov et 
al. m who studied the Jeans instability of a relativistic SF in a static background. The case of a 
complex self-interacting SF representing BECDM was considered by Chavanis m in the context 
of Newtonian cosmology. His study is based on the GPP equations. Harko m and Chavanis m 
developed an approximate relativistic BEC cosmology and found that the perturbations grow 
faster in BECDM as compared to ACDM. Recently, Li et al. da developed an exact relativistic 
cosmology for a complex self-interacting SF/BEC based on the KGE equations. They studied 
the evolution of the homogeneous background and showed that the Universe undergoes three 
successive phases: a stiff matter era, followed by a radiation era due to the SF (that exists only 
for self-interacting SFs), and finally a matter era similar to CDM. 

Instead of working directly in terms of a SF, we can adopt a fluid approach and work with 
hydrodynamic equations. In the case of the Schrodinger equation, this hydrodynamic approach 
was introduced by Madelung [18] . He showed that the Schrodinger equation is equivalent to the 
Euler equations for an irrotational fluid with an additional quantum potential arising from the 
finite value of Ti. This hydrodynamic representation has been used by Bohmer and Harko m 
and by Chavanis pa EH El] among others in the case of BECDM and in the case of BEC stars 
|22j. This hydrodynamic approach has been generalized by Suarez and Matos [12, [23j in the 
context of the KG equation. They used it to study the formation of structures in the Universe, 
assuming that DM is in the form of a fundamental SF with a Ay? 4 potential. 

In the works P2E3], the SF is taken to be real and the gravitational potential is introduced 
by hand in the KG equation, and assumed to be determined by the classical Poisson equation 
where the source is the rest-mass density p. This leads to the KGP equations. However, this 
treatment is not self-consistent since it combines relativistic and nonrelativistic equations. In 
these Proceedings, we derive the hydrodynamic representation of a complex SF coupled to 
gravity through the Einstein equations in the weak field approximation. In this way, we develop 
a self-consistent relativistic treatment. Throughout this work, we use the conformal Newtonian 
gauge which takes into account metric perturbations up to first order. We consider only scalar 
perturbations. This is sufficient if we are interested in calculating observational consequences 
of the SF dynamics in the linear regime. We compare our results with those obtained from the 
heuristic KGP equations. In these Proceedings, we develop the basic formalism and provide 
the main equations for an arbitrary self-interaction potential of the form U(|y?| 2 ). Details and 
applications of our formalism can be found in our research papers [ 231 125] . 


2. The conformal Newtonian gauge 

We consider the KGE equations in the weak field limit < F/c 2 -C 1. The equations that we 
derive are valid at the order (){<&/(?). We work with the conformal Newtonian gauge which is 
a perturbed form of the Friedmann-Lemaitre-Robertson-Walker (FLRW) line element [26]. We 
consider the simplest form of the Newtonian gauge, only taking into account scalar perturbations 
which are the ones that contribute to the formation of structures in cosmology. We neglect 
anisotropic stresses. We assume that the Universe is flat in agreement with the observations of 
the cosmic microwave background (CMB). Under these conditions, our line element is given by 


ds z = c 2 ( 1 + 2-^- ) dt z — a(t ) 


$ 

1 — 2— ) 5ijdx l dx 3 , 


(1) 


where <f>/c 2 <C 1. In this metric, <h(x, t) represents the gravitational potential of classical 
Newtonian gravity and a(t) is the (dimensionless) scale factor. 


3. The Lagrangian of the scalar field 

We assume that DM can be described by a complex SF which is a continuous function of space 
and time defined at each point by p(x = p(x,y, z,t). The action of the relativistic SF is 

S v = J ( 2 ) 

where = C^p, p*, B^p, B^p*) is the Lagrangian density and g = det(g^ u ) is the determinant 
of the metric tensor. We adopt the following generic Lagrangian density 

M 2 -r(M 2 ) (3) 

which is written for a metric signature The second (quadratic) term is the rest- 

mass term. The potential E(|t/?| 2 ) takes into account the self-interaction of the SF. In certain 
applications, it is relevant to consider a quartic potential of the form 

V{\V\ 2 ) = §^M 4 , (4) 

where A is the self-interaction coupling strength. If the SF describes a BEC at T = 0, the 
quartic potential ([5]) models the two-particle self-interactionQ In that case, m represents the 
mass of the bosons and A is a constant related to the s-wave scattering length of the bosons a s 
(measured in the center of mass frame) by A = 4:7ra s h 2 /m. This corresponds to the first Born 
approximation. The potential can then be written as 

Tr/I .n. 27 t a s m. .a 

V(\p\ 2 ) = ( 5 ) 


4. The energy-momentum tensor 

The energy-momentum tensor of the SF is 

BC, - 

= ~ ( 6 ) 

For the generic Lagrangian Q, it takes the form 


= -(d^p*d v p + d u p*dfj,p) - g^ 


m?c 2 


^g pa d p p*d a p-— T \p\ — V{\p\ ) 


( 7 ) 


By analogy with the energy-momentum tensor of a perfect fluid, the energy density and the 
pressure tensor of the SF are defined by e = Tg and P- = — T/. The conservation of the 
energy-momentum tensor, which results from the Noether theorem, writes D U T IIU = 0. 

1 It is a good approximation to ignore higher-order interactions when the boson gas is dilute, i.e., when the 
particle self-interaction range is much smaller than the mean interparticle distance. 








5. The Klein-Gordon equation 

The equation of motion for the SF can be obtained from the principle of least action. Imposing 
5Sip = 0 for arbitrary variations 5p and dip*, we obtain the Euler-Lagrange equation 


D„ 


dC u 


d{d^y 


dC, 

dip* 


= 0, 


( 8 ) 


where D is the covariant derivative. For the generic Lagrangian this leads to the KG 
equation 


m 2 c 2 


H- ip + 2K(|(/?| ),y>* — 0, 


where □ is the d’Alembertian operator 


and 


D = D^g^dv) = y==<9 /i (\/—5 g^ v d u ) 




(9) 

( 10 ) 

( 11 ) 


Computing the d’Alembertian (1101) with the Newtonian gauge, we obtain the KG equation 


1 d 2 p^3H dip 1 


c 2 dt 2 


dt 


1 + 


4<h 


A p— 


4 d$ dp 
c 4 dt dt 


+ 1 + 


24>\ m 2 c 2 
c 2 ) h 2 


<p+2(l + 2^jv,v* = 0, ( 12 ) 


where H = a/a is the Hubble constant. Using the expression (J7| of the energy-momentum 
tensor, the energy density and the pressure are given by 


f — T u — 
e — i 0 — 


1 

2c 2 


1 - 


24>\ 

I?) 


dp 


dt 


+ 


1 

2^ 


1 + 


24* 


,2 _2 


\vp\ z +^\p\ 2 +n\pn 


(13) 


p = -i(r 1 1 + T 2 + T 3 3 ) 


J_ (, _ 2£\ dp 2 

2 c 2 \ c 2 ) dt 


1 

6a 2 



m 2 c 2 

2h 2 


M 2 


V(\p\ 2 ). 

(14) 


6. The Einstein equations 

The Einstein equations write 


RflV 


8 

- T 

c 4 1 V’ v ' > 


(15) 


where R pu is the Ricci tensor, R is the Ricci scalar, and G is Newton’s gravitational constant 
m- The conservation of the energy-momentum tensor is automatically included in the Einstein 
equations. The time-time component of the Einstein equations is 


e>° lp_ 87rG o 
-On — -a — -rb ■ 


(16) 


With the Newtonian gauge, it can be written as 


A 4* e 3 H 2 3 H ( S4> TT ^ 

-=-1-(-h P4> 

47rGa 2 c 2 8nG 4 ttGc 2 v dt 


( 17 ) 
























Using the expression (11311 of the energy density e, which represents the time-time component of 
the energy-momentum tensor, we get 


A4> 

\txGo? 


1 

2? 


1 - 


2$ 


dip 


dt 


+ 


2 a 2 c 2 


1 + 


2<F\ 


m 


1 .. l2 , 3 H 2 

+-»v{ M 2 )- —+ 


c 2 J 

3 H 




87 tG AitGc 2 


<94> 

aF + ff,|> 


Eqs. d and (fT 8 l) form the KGE equations. 


(18) 


7. Spatially homogeneous scalar field 

For a spatially homogeneous SF with ipb(x,t) = <pb{t) and t) = 0, the KG equation (fl2l) 
reduces to 


1 d 2 ip b 3 H_dp>b_ 
c 2 dt 2 c 2 dt 


m 2 c 2 

IF" 


+ 2U j¥ ,* 


= 0. 


(19) 


In that case, the energy-momentum tensor is diagonal and isotropic, Tj^ = diag(e&, —Pb, — Pb , —Pb)- 
The energy density e&(t) and the pressure P;,(t) are given by 


€b = 


2 c 2 


difb 


dt 


+ 


m 2 c 2 
2h 2 


\V>\ 2 + V{\<Pb\ 2 ), 


Pb = 


2 c 2 


difb 


dt 


m 2 c 2 

2h 2 




( 20 ) 


From these equations, we obtain the continuity equation 

+ A) = 0 (21) 

which is one of the Friedmann equations m- The other Friedmann equation is obtained from 
the Einstein equation (fTKl) that reduces to 


H 2 


8 v tG 
"3c 2 


-£b- 


( 22 ) 


This relation shows that the term —3H 2 /8ttG = —eb in the Einstein equation (11711 plays the role 
of a neutralizing background!! From Eqs. (12111 and (|22|1 . we easily obtain 


a 

a 


4t tG 
"3c 2 " 


( e 6 + 3-Pfr). 


(23) 


8. The Gross-Pitaevskii-Einstein equations 

The KG equation without self-interaction can be viewed as a relativistic generalization of the 
Schrodinger equation. Similarly, the KG equation with a self-interaction can be viewed as a 
relativistic generalization of the GP equation. In order to recover the Schrodinger and GP 
equations in the nonrelativistic limit c —>■ + 00 , we make the transformation [I]: 

<p(x, t ) = -e- imc2t/7 V(z, t). (24) 

m 

2 In a static universe, the source of the gravitational potential is e/c 2 . In an expanding universe, the source of 
the gravitational potential is (e — £b)/c 2 , so that $ = 0 when e = tb- When we work in the comoving frame, the 
expansion of the universe amounts to subtracting a neutralizing background — et/c 2 to the density e/c 2 [25j, as 
in the Jellium model of plasma physics. This is the correct way to solve the Jeans swindle [29]. 






















The prefactor Ti/m is justified in Appendix A Mathematically, we can always make this change 
of variables. However, we emphasize that it is only in the nonrelativistic limit c —>• +oo that ip 
has the interpretation of a wave function, and that \ip\ 2 = p has the interpretation of a rest-mass 
density. In the relativistic regime, ip and p = \ip\ 2 do not have a clear physical interpretation. 
We will call them “pseudo wave function” and “pseudo rest-mass density”. Nevertheless, it is 
perfectly legitimate to work with these variables and, as we shall see, the equations written in 
terms of these quantities take relatively simple forms that generalize naturally the nonrelativistic 


ones. 


Substituting Eq. 


into Eqs. m and (j!8D , we obtain 


dip _ h 2 d 2 ip 
dt 2me 2 dt 2 


- -H 


n 2 dip 
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h 2 


- 1 + —^\m 


2 ma? 
2 


4$' 

1 H—— ) Aip — m&ip 


2 me 2 dt 

24* \ dV , 3 , , 2fi 2 <94> ( dip ime 


iip + -iUHip + 


d\ip\ 2 2 me 


4 dt 


dt h w °’ 


(25) 


A4> 


AirGa 2 


= n-* 


+ 


n 2 


2m? c 4 


+V(h/f)- 


me* 


2$\ 

1 C^" ) Im 


2<f \ 

1 *^) 

dip 


dip 


dt 


+ 


n 2 


2 a?m 2 c 2 


A 3 H 2 3 H 

dt ) 8 itG 47tGc 2 V dt 


i + ^’)lAl 2 

<94* 


+ H4* 


(26) 


Eq. (1251) can be interpreted as a generalized Schrodinger equation (in the absence of self¬ 
interaction) or as a generalized GP equation (in the presence of self-interaction). It is coupled 
to the Einstein equation ([26D . The energy density and the pressure can be written as 


n 2 


e = 


2 m 2 c 2 


24* 

1-7T" 


24* 


^ 2 _ 1U _ A ) Im 

m \ c 2 / 


dt 


dt 


24* 


+2A? I 1 + ^ ™+ rT - w + -pw 2 + V(m, 


1 


(27) 


P = 


h 2 


2m? c 2 


1 - 


24* 


dip 2 h 
~dt 


m v c z / \ at 


w (i + ?) 'Ai 2 + r 2 (i - f) ivf - - V(\p). 


(28) 


Eqs. (1251) and (j26[) form the GPE equations. In the nonrelativistic limit c —>• +oo, they reduce 
to the GPP equations [15]: 


.. dip 3 fi 2 dV 

hH* = + m** + 


(29) 


A4* 

47 t Ga? 

For the quartic potential we have 


3 H 2 
8nG' 


(30) 


vm 2 ) = 


2Tra.fi 2 




m 4 


(31) 

































9. The hydrodynamic representation 

Important characteristics of the system are revealed by reformulating the KGE equations in the 
form of hydrodynamic equations. This can be done at the level of the GPE equations (j25H - (126(1 
via the Madelung transformation [18]. To that purpose, we write the pseudo wave function ijj 
as 

^(x,t) = ^J)e iS ^ h , (32) 

where p = |^| 2 plays the role of a pseudo rest-mass density and S plays the role of a pseudo 
action. Following Madelung, we also define a pseudo velocity field as 


v(x, t) 


VS 
ma ’ 


(33) 


where the scale factor a has been introduced in order to take into account the expansion of the 
Universe. We note that this velocity field is irrotational. 

Substituting Eqs. (l32l) - (|33l) into the GPE equations (l25l) - (l26l) . and separating real and 
imaginary parts, we obtain the system of hydrodynamic equations 




J_d_( 9S\ 3H_pdS_ _4/^d£ / 2 

me 2 dt v dt) + me 2 dt me 4 dt \ 


dS\ 

~dt) 


4$ - 


ac* 


(34) 


dS (VS) 2 
dt 2 ma 2 



9 2 Vp 


2 me 2 y/p 
mh(p ) + 


+ 1 + 


44>\ n 2 Ay/p 


24> 


1 


2 me 2 V dt J 


c 2 / 2 ma 2 
dS x2 


•V- 


3 H - 


VP 

4 d4>\ 


me 2 a 2 


(VS) 5 


n 2 dp 


c 2 dt J Amc 2 p dt ’ 


(35) 


dv 1 -> h 2 

- + Hv + - (v . V ) v = -—— 


V 


~dW 

Vp 


+ 


h 2 


2 m 2 a? 


V 


c 2 J Vp 


1 


--V$- —Vp- -AjV(/i$) - ^V($u 2 ) + o 

a pa ac - ac z 2am 2 c 2 


V 


dSV 

dt 7 / 


3ft 2 ^ , 

4am 2 c 2 \p dt 


h 2 


am 2 c 4: 


V [ 

dt p dt 


A4> 

47rGa 2 


= 1 - 




P + 


n 2 


+ 


2m 2 c 4 
ft 2 


24>\ 


2 a 2 m 2 c 2 


1 + 


2<h 


i® 2 , p® 2 
4 p \dt) n 2 \dt) 

1 


+\v(p) - 

cr mc z 


2$\ dS 3H 2 
~^r) P ~dt ~ 8 ttG + 4v tGc 2 


(Vp) 2 + ^(VS) 2 
3 H 


4p' ' r ' ' ?i 2 


d<k 


(36) 


(37) 


where /i(p) = V'(p) is a pseudo enthalpy and p(p) is a pseudo pressure defined by the relation 
h!(p) = p'{p)/p pnj- It is explicitly given by p(p) = ph(p) — J h(p ) dp, i.e., 


p(p) = pV'(p) -V(p). 


( 38 ) 



































The pseudo velocity of sound is c 2 = p'(p) = pV"{p). For the quartic potential ©, we have 

2-k a s h 2 2 47ra s ?i 2 2na s h 2 2 9 4Tra s h 2 

= -3 —P , Kp) = -3 —P, P{p) = -3 —P , c s = -3 —P- 

m J m 6 m J m 6 

The pseudo pressure is given by a polytropic equation of state of index 7 = 2 which is quadratic. 
We note that this equation of state coincides with the equation of state of a nonrelativistic self¬ 
interacting BEC [30]. This coincidence is not obvious because Eqs. (|MD - (13T1) are valid in the 
relativistic regime. The interpretation of this equation of state is, however, not direct because 
p and p are a pseudo density and a pseudo pressure that coincide with the real density and the 
real pressure of a BEC only in the nonrelativistic limit c —>• + 00 . 

The energy density and the pressure can be written in terms of hydrodynamic variables as 




h 2 f 24>\ 

2 m 2 c 2 V c 2 ) 
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p fdS V 


h 2 \dt J 
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(Vp) 2 + ^(V5) 
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pc +-pc +V(p), 


(40) 


P = 


n 2 


2 m 2 c 2 


2$ 
1-77- 


J_ f^P \ 2 , P_ 


4 p\dt J + h 2 \dt J 


6 a 2 




, 24* 

H— 2 " 

c z 


± ( V P f + £(VSf 


, 2$\ p dS 1 / 2$\ 2 1 2 rt N 

)^^+o (1- —)pc -7PC -E(p). 


m 




(41) 


We note that, in general, the pressure P defined by Eq. m is different from the pressure p 
defined by Eq. (13811 FI However, they coincide for a homogeneous SF in the regime where the SF 
oscillations are faster than the Hubble expansion (24, 25]. 

The hydrodynamic equations (|34D - (137|1 have a clear physical interpretation. Eq. (1341) . 
corresponding to the imaginary part of the GPE equations, is the continuity equation. We note 
that / p d 3 x is not conserved in the relativistic regime. However, Eq. (1341) is consistent with the 
conservation of the charge of a SF (see section flOl and Appendix B ). Eq. ([35]) . corresponding 
to the real part of the GPE equations, is the Bernoulli or Hamilton-Jacobi equation. Eq. (1361) . 
obtained by taking the gradient of Eq. (13511 . is the momentum equation. Eq. (13711 is the 
Einstein equation. We stress that the hydrodynamic equations (l341) - (l3Tl) are equivalent to the 
GPE equations (f25|) - ([26j) which are themselves equivalent to the KGE equations (fT2|) and (fTKl) . 
In the nonrelativistic limit c —>• + 00 , we recover the quantum Euler-Poisson equations [T5 j: 


% + 3Hp+-V-(pv) = 0 , 
at a 


(42) 


dS_ {VS) 2 
dt 2 ma 2 


ft 2 A Jp 

2 ma 2 y/p 


— m4> — mh(p), 


d y TT 1 
— + Hv + -{v ■ V)v = 
dt a 


h 2 

2 m 2 a 3 



— -V4>-— Vp, 

a pa 


(43) 

(44) 


A4* _ 3 H 2 

AnGa 2 ^ 87 xG 

These equations can also be written in terms of the density contrast 5 = (p — Pb)/Pb [IS]. 

3 In our formalism, p represents the pressure arising from the self-interaction of the bosons (scattering) while P 
is the pressure of the SF defined by analogy with the pressure of an ideal fluid whose energy momentum tensor 
writes TA = (P + t)u^u v j(? - Pg 


























10. Cosmological evolution of a spatially homogeneous scalar field 

We consider the evolution of a universe induced solely by a spatially homogeneous SF. In the 
comoving frame, we have p(x,t ) = pb(t), v b (x,t ) = 0, 3>&(x, £) = 0, and S b {x,t) = Sb(t). 
We introduce the notation E(t) = —dSb/dt which can be considered as the time-dependent 
energy of the spatially homogeneous SF in the comoving frame. The pseudo wave function 
of the SF is ^(x,£) = iJjb(t) = \fpbif)e~^^ fl ^ E ( t ) dt . Using Eq. (ITU) . we have tpb(x,t) = 
Mt) = A^t) e -m)[mch + jE^dt] so the total energy of the SF, including its rest mass, is 
Etat{t) =E(t) + me 2 . 

For a spatially homogeneous SF, the hydrodynamic equations (1M1) - (07D reduce to 


dpb 

dt 


+ 3 Hpb — 


1 d 
me 2 dt 


(. PbE ) 


3 Hp b 
me 2 


E, 


(46) 
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2 me 2 



E = 


* 2 d2 VPb 
n -dE~ 

2 me 2 yfpb 


+ mh(p b ) + 


3 Hh 2 dpb 
Amc 2 p b dt 


(47) 


3 H 2 
8ttG 


ft 2 

^ b + 2 m 2 c 4 


1 (dpbV ,fE F 2 

4 pb V dt ) h 2 


1 Jr , . E 
+ (W “I-2^ 6 ' 

mc z 


(48) 


In terms of the total energy E tot (t) = E(t ) + me 2 , the equation of continuity (HS1) becomes 


1 dp b 3da 1 dE tot 

Pb dt a dt E tot dt 

It can be rewritten as a conservation law: 

^ (Et 0 t Pba 3 ) = 0 . 

Therefore, the total energy is exactly given by 

Etot _ Qm 
me 2 pbd 3 ’ 


(49) 


(50) 


(51) 


where Q is a constant. This conservation law was found by Gu and Hwang m directly from 
the KG equation. It can be shown that Q = J J° d 3 x represents the conserved charge density of 
the complex SF (see |Appcndix B| ). The energy density and the pressure of a homogeneous SF 
are 


£b 


h 2 1 

8 m 2 c 2 pb 


(dp b \ 2 Pb 
\ dt ) 



E \ 

2me 2 ) ~*~ 


p b c 2 + V(p b ), 


(52) 


Pb 


n ' 2 1 (^Pb \ 2 , Pb E ( x , E 

8 m 2 c 2 pb V dt J m \ 2 me 2 


-V(pb). 


(53) 


Equations m, ssi) and m determine the complete evolution of a universe induced by a 
spatially homogeneous SF. Working directly on the homogeneous KG equation with a quartic 
self-interaction potential, Li et al. m have shown that a universe filled with a relativistic 
complex SF first undergoes an intrinsic stiff matter era, followed by a radiation era due to its 
self-interaction, before finally entering in the matter era. The stiff matter era occurs when the 
























SF oscillations are slower than the Hubble expansion while the radiation and matter eras occur 
when the SF oscillations are faster than the Hubble expansion. These different regimes can be 
recovered from the hydrodynamic equations (14611 - (I48|) }24j . 

In the nonrelativistic limit c —» +oo, Eqs. (146D - ()48(i reduce to 

^- + 3Hp b = 0, E = mh(p b ), = Pb- ( 54 ) 

We find that p b oc a” 3 , a oc i 2 / 3 and p b = l/( 67 rGi 2 ) (Einstein-de Sitter solution), so 
the homogeneous SF/BEC behaves as CDM. For the quartic potential ©, we have E(t) = 
Airash 2 p b /m 2 = 2a s h 2 /3Gm 2 t 2 and S b (t) = 2a s h 2 /3Gm 2 t + C. 


11. Generalized Klein-Gordon-Poisson equations 

In this section, we consider a simplified model in which we introduce the gravitational potential 
<h(x, t) in the ordinary KG equation by hand, as an external potential, and assume that this 
potential is produced by the SF itself via a generalized Poisson equation in which the source is 
the energy density e. This leads to the generalized KGP equations. We then show that these 
equations can be rigorously justified from the KGE equations in the limit <f>/c 2 —>• 0. However, 
this simplified model is not sufficient to study the evolution of the perturbations in the linear 
relativistic regime since it precisely neglects terms of order <h/c 2 . 

We consider the FLRW metric that describes an isotropic and homogeneous expanding 
background. The line element in the comoving frame is 

ds 2 = g lw dx ,l dx lJ = c 2 dt 2 — a{t) 2 5ijdx l dxE (55) 


For this metric, the d’Alembertian operator ( 1101 ) writes 


Id 2 3 H d 1 
c 2 dt 2 c 2 dt a 2 


(56) 


In order to take the self-gravity of the SF into account, we introduce a Lagrangian of interaction 
that couples the gravitational potential to the scalar field ip(x,t) according to 


Ant = 


(57) 


The total Lagrangian of the system (SF + gravity) is given by £ = + C ni t. The equation 

of motion resulting from the stationarity of the total action S = S v + 5; n t, obtained by writing 
<55 = 0, is the KG equation 

y-Y-f 2 ^>2 O-vvi 2 

4 ¥ + 2K(M 2 )j¥3* 4 = 0, (58) 


where the d’Alembertian operator is given by Eq. (1561) and the gravitational potential 4>(x, t) 
acts here as an external potential. The energy density and the pressure, defined from the 
diagonal part of the energy-momentum tensor ( 0 , are given by 
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Eq. (1581) is the ordinary KG equation for a SF in an external potential d*(x, t) in an expanding 
background. We now state that d>(x, f) is actually the gravitational potential produced by the 
SF itself. We phenomenologically assume that the gravitational potential is determined by a 
generalized Poisson equation of the form 

Ad* 1 

= < 61 > 

in which the source of the gravitational potential is the energy density e of the SF (more precisely, 
its deviation from the homogeneous background density e&(i)). Using Eq. (1591) for the energy 
density of a SF, and recalling the Friedmann equation (I22p . the generalized Poisson equation 
can be written as 


Ad* 

AirGa 2 
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2c 1 


d(p 
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(62) 


Eqs. (1581) and (1621) form the generalized KGP equations. They have been introduced in an ad 
hoc manner but they can be rigorously justified from the KGE equations (1121) and (1181) in the 
limit d>/c 2 —t- 0 (which, of course, is different from the nonrelativistic limit c —> +oo). We see 
that the gravitational potential d> appears in the KG equation (1581) due to the cancelation of 
c 2 in the product d*/c 2 x c 2 in Eq. (1121) . Therefore, we do not have to introduce d> by hand: 
the generalized KGP equations can be obtained from the KGE equations by simply neglecting 
terms of order d>/c 2 in these equations. Similarly, the equations related to the generalized KGP 
equations can be obtained from the ones related to the KGE equations by neglecting terms of 
order d>/c 2 . For example, the generalized GPP equations write 
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The energy density and the pressure are given by 
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The corresponding hydrodynamic equations write 
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The energy density and the pressure can be written in terms of hydrodynamic variables as 
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-V(p). (72) 


This model correctly describes the homogeneous background for which <f> = 0 but it is not 
sufficient to describe the evolution of the perturbations in the linear regime because we must 
precisely take into account the terms of order <b/c 2 in this regime (except, of course, in the 
nonrelativistic limit c —>• +oo). Therefore, the use of the KGE equations is mandatory to study 
the evolution of the perturbations in the relativistic regime. 

Remark: We could also assume that the gravitational potential is determined by a Poisson 
equation of the form 


Ad* 

AirGa 2 ^ 


(73) 


in which the source of the gravitational potential is the pseudo rest-mass density p = \ip\ 2 of 
the SF. Eqs. (1581) and (1731) form the KGP equations. This approximation has been considered 
hr 52j. However, there is an inconsistency in coupling the relativistic KG equation (1581) to the 
classical Poisson equation (1731) . 


12. Conclusion 

We have developed a formalism based on a relativistic SF described by the KGE equations 
in the weak held limit. We have transformed these equations into equivalent hydrodynamic 
equations. These equations are arguably more tractable than the KGE equations themselves. In 
the nonrelativistic limit, they reduce to the hydrodynamic equations directly obtained from the 
GPP equations m- Therefore, our formalism clarifies the connection between the relativistic 
and nonrelativistic treatments. We note that, in the relativistic regime, the hydrodynamic 
variables ip, p, S, v, p,... that we have introduced do not have a direct physical interpretation. 
It is only in the nonrelativistic limit c —>• +oo that they coincide with the wave function, rest- 
mass density, action, velocity, and pressure. However, these variables are perfectly well-defined 
mathematically from the SF p in any regime, and they are totally legitimate. Furthermore, in 
terms of these variables, the relativistic hydrodynamic equations take a relatively simple form 
that provides a natural generalization of the nonrelativistic hydrodynamic equations. 

The complete study of these relativistic hydrodynamic equations is of considerable interest but 
it is, of course, of great complexity. In our research papers [MIES], we have started their study in 
simple cases. We have checked that the hydrodynamic equations of the SFDM model reproduce 
the evolution of the homogeneous background obtained previously by Li et al. m directly 
from the KGE equations: a stiff matter era, followed by a radiation era (for a self-interacting 
SF), and a matter era. We have also started to study the evolution of the perturbations in 
the linear regime in a static and in an expanding universe. We have shown analytically and 
numerically that perturbations whose wavelength is below the Jeans length oscillate in time 
while perturbations whose wavelength is above the Jeans length grow linearly with the scale 
























factor as in the CDM model. The growth of perturbations in the SF model is substantially 
faster than in the CDM model. We have also shown that general relativity attenuates or even 
prevents the growth of perturbations at very large scales, close to the horizon (Hubble length). 

For physically relevant wavelengths, the nonrelativistic limit of our formalism is sufficient 
to describe the evolution of the perturbations in the matter era. However, even if relativistic 
corrections are weak in the matter era, we may wonder whether their effect could be detected 
in an era of precision cosmology. In particular, it would be interesting to see if one can observe 
differences between the KGE equations considered in these Proceedings and the heuristic KGP 
equations studied in the past, in which gravity is introduced by hand in the KG equations (see 
section HU). 

In future works, it will be important to consider the nonlinear regime where structure 
formation actually occurs. In general, this problem must be addressed numerically. The 
hydrodynamic equations derived in these Proceedings may be very helpful because they may be 
easier to solve than the KGE equations. Therefore, numerical simulations using fluid dynamics 
should be developed in the future. As a first step, relativistic effects could be neglected and 
the nonrelativistic equations of Ref. m could be used. These equations are similar to the 
hydrodynamic equations of CDM except that they include a quantum potential (Heisenberg) 
and a pressure term (scattering) that avoid singularities at small scales |15J . In this respect, it 
may be recalled that the SP equations were introduced early by Widrow and Kaiser [32] as a 
procedure of small-scale regularization (a sort of mathematical trick) to prevent singularities in 
collisionless simulations of classical particles. In their approach, h is not the Planck constant, 
but rather an ajustable parameter that controls the spatial resolution. Their procedure may 
find a physical justification if DM is made of self-gravitating BECs [33J. 

Our relativistic formalism may have applications for other self-gravitating systems described 
by SFs or BECs besides DM. We can mention, for example, the case of boson stars [3EUZ] and 
the case of microscopic quantum black holes made of BECs of gravitons stuck at a quantum 
critical point [3H [351 HO]- It has also been proposed [22] that, because of their superfluid core, 
neutron stars could be BEC stars. Indeed, the neutrons (fermions) could form Cooper pairs and 
behave as bosons of mass 2 m n . This idea may solve certain issues regarding the maximum mass 
of neutron stars. Finally, we may mention analog models of gravity in which BECs described by 
the GP equation or by the KG equation are used to simulate results of classical and quantum 
held theory in curved spacetime [36]. 
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Appendix A. The value of A 

The GP equation is obtained from the KG equation by means of the transformation 

tp = Ae- irm?t / n ip. (A.l) 

The constant A can be computed as follows. Substituting Eq. (IA.1I) into Eq. (fT3l) . we find that 
the energy density of the SF is given by 
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find 


is the rest-mass density. Since e ~ pc 2 in the nonrelativistic limit c —>• +oo, we 
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Appendix B. Some comments about the Klein-Gordon equation 

In this Appendix, we recall the difficulties associated with the interpretation of the KG equation. 
We also clarify the relation between the hydrodynamic representation of the KG equation and 
the conservation of the charge. 

The fundamental equation of nonrelativistic quantum mechanics is the Schrodinger [3] 
equation 




(B.l) 


If we write the wave function under the form ip(x,t) = \Jp(x, t)e l ‘ an d define the density 
p and the current J by 


P = lv>l 2 , 


(B.2) 


J = P™ = 2 ^*^ - W*), (B.3) 

where we have used S = (h/2i) we obtain the continuity equation 

^ + V-J = 0. (B.4) 

This equation shows that the integral of the density J pd?x is conserved. Furthermore, by 
definition, the density is positive: p(x, t ) > 0. Therefore, p(x, t ) can be interpreted as a 
probability density. 

In the Madelung hydrodynamic representation of the Schrodinger equation, we define the 
density p and the velocity v by 


p=m , 

These variables satisfy the continuity equation 
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Since J = pv according to Eqs. ()B.3D and (1B.5I) . we immediately see the equivalence between 
Eqs. (IB. 41) and (IB.61) . 

We now consider the KG mm equation 
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which was initially proposed as a relativistic extension of the Schrodinger equation^ If we 
write the SF under the form <p(x, t ) = \JR(x , t)e ia< ^ x,t ^ n and introduce the quadricurrent 


4 The KG equation was actually discovered by Schrodinger before he found the equation that now bears his name 
DO- The KG equation was also obtained by Fock [37], de Donder and van den Dungen [38], and Kudar [39]. 











= —Rd^cr/m = —(h/2im)(p*d fJ, p — pd^p*), we obtain the continuity equation d^J^ = 0. 
The quadricurrent can be written as J^ = (J°, J) with 
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J = R- = —(p*Vp-pVp*). (B.9) 

m 2 im 

If we introduce the KG density pkg = J° /c, we can write the continuity equation as 

%°+V-J = 0. (B.10) 

dt y ’ 

This equation shows that the integral of the KG density J pkqc1 3 x is conserved. However, 
the KG density pKc(x,t) is not definite positive so it cannot be interpreted as a density 
probability. Another difficulty with the KG equation is that it allows negative kinetic energies 
as solution. Indeed, decomposing Eq. (IB. 71) into plane waves, we obtain two solutions of the 
form p±(x,t) = A±e l (P' x ~ E± dl n with E± = ±\/p 2 c 2 + m 2 c 4 . 

The original difficulties encountered with the KG equation had interesting historical 
developments. Dirac [HI] proposed another relativistic extension of the Schrodinger equation. 
The Dirac equation describes spin-1/2 massive particles (fermions) such as electrons and quarks. 
In Dirac’s theory, the probability density is positive but negative energies are allowed. Dirac 
solved the problem of negative energies via the “hole” theory SU- This leads to the concept of 
antiparticles that are related to negative energy eigenstates. Antimatter was unsuspected before 
Dirac’s work. The first antiparticle, the positron, was experimentally discovered by Anderson 
[42] in 1932. On the other hand, Pauli and Weisskopf [43] proposed a new interpretation of the 
KG equation. They interpreted the KG density pkg = J°/c as a charge density which can be 
of arbitrary sign. If we define the charge density and the charge current by p e = meJ°/ch 2 and 
J e = meJ/h 2 , where e is an elementary charge, we can rewrite the continuity equation (IB.101) 
as 

^ + V-J e = 0. (B.ll) 

This equation expresses the conservation of the total charge Q e = J p e d 3 x. Although the 
KG equation is not a successful relativistic generalization of the Schrodinger equation (initially 
introduced to describe the energy spectrum of the electron), this equation was resurrected in the 
context of quantum held theory where it was shown to describe spin-0 particles (bosons) such 
as 7r-mesons, pions, or the Higgs boson. In this interpretation, since p^ = ±(e\E±\/Ti 2 c 2 )\p±\ 2 , 
p+ specifies particles with charge +e and energy +E while p- specifies antiparticles with the 
same mass but with charge — e and energy —E. For a real SF the charge is zero since p e = 0 
and J e = 0 according to Eqs. (IB. 81) and (IB.91) . 

In the hydrodynamic representation of the KG equation, we write 

p(x,t) = -Jp(x,ty\ s{ - 3 ^- mc2t V n (B.12) 
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These variables satisfy the continuity equation 
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The density p(x, t) is positive by definition but it is not conserved. Therefore, p(x, t ) cannot be 
interpreted as a density probability. However, the continuity equation (IB.141) can be rewritten 
as 


d_ f E tot \ 
dt v me 2 ) 


+ V • (pv) = 0, 
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where we have defined E(x,t ) = —dS/dt and E tot (x,t) = me 2 + E(x,t). This equation implies 
the conservation of the integral / pE tot d 3 x. To show that this integral corresponds to the total 
charge of the SF, we note that R = {Ti 2 /m 2 )p and a = S — mc 2 t. Therefore, 
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These relations first establish the equivalence between Eqs. (IB. 111) . (IB.14j) . and (IB.15p . 
Furthermore, they show that the total charge of the SF can be written as 


Q e = e [ ——^ d 3 x (B.18) 

J m mc z 

in agreement with Eq. (IB.151) . Eq. (IB. 181) is consistent with the expression (1511) of the charge 
of a spatially homogeneous SF in an expanding universe. 
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